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Abstract. Adding quenched disorder to the one-dimensional asymmetric exclusion 
process is known to always induce phase separation. To test the robustness of this 
result, we introduce two modifications of the process that allow particles to bypass 
defect sites. In the first case, particles are allowed to jump I sites ahead with the 
probability pi ^ where a > 1. By using Monte Carlo simulations and the 

mean-field approach, we show that phase coexistence may be absent up to enormously 
large system sizes, e.g. InL ^ 50, but is present in the thermodynamic limit, as in the 
short-range case. In the second case, we consider the exclusion process on a quadratic 
lattice with symmetric and totally asymmetric hopping perpendicular to and along the 
direction of driving, respectively. We show that in an anisotropic limit of this model a 
regime may be found where phase coexistence is absent. 

Keywords: driven diffusive systems (theory), disordered systems (theory), stationary 
states 
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1. Introduction 

Disorder usually lias a strong impact on critical phenomena, a hallmark example being 
the influence of spatially uncorrelated quenched disorder on the second-order phase 
transition, the relevance of which is predicted by the Harris criterion [1]. Such a universal 
result is, however, generally lacking in systems far from equilibrium. There are few 
exceptions, and particularly well-studied are driven diffusive systems, in particular the 
asymmetric simple exclusion process (ASEP), for which even a version of the Harris 
criterion can be established [2]. Being exactly solvable [3], ASEP has attracted a lot 
of attention due to its connection with quantum spin chains nonequilibrium phase 
transitions [5] and wide range of applications like surface growth [6], traffic [7] and 
biological transport [8]. 

ASEP is a simple model describing classical particles hopping in the preferred 
direction on a discrete lattice and interacting only by the exclusion principle that forbids 
two particles from occupying the same lattice site (for reviews, see e.g. [Hllin])- Disorder 
in ASEP is introduced by assigning random (but quenched) hopping rates either to sites 
or particles. In the simplest case of a single defect particle hopping at rate r < 1, the 
exact solution was found [HI [T2], with the conclusion that the defect induces phase 
coexistence through the creation of shocks for all particle densities p < 1 — r. A 
similar conclusion was obtained for a model with an arbitrary number of slow particles 
solved simultaneously by Evans [13] and Ferrari and Krug [H], who observed that the 
appearance of shocks is determined solely by the slowest particle in the system. On the 
contrary, ASEP with a single defect site is still an open problem. Although early studies 
of this problem in various contexts [151 UHl [IZl [HI [I9] provided strong arguments that 
a single defect always induces phase separation (i.e. for all r < 1), several later studies 
suggested the existence of a threshold value rc < 1 above which the shocks disappear 
[201 [2T| . On the other hand, no such "shock- free" scenario is possible for the full site- 
wise disorder, as pointed out by Tripathy and Barma [221 123] in the particular case of a 
binomial distribution. The reason behind this is essentially a geometric one, in a sense 
that disorder creates a long stretch ( "bottleneck" ) of slow sites that limits the current 
and ensures the shocks to appear as soon as this limiting value is reached. 

The fact that phase coexistence appears for any distribution of disorder [21] raises 
the following question: do such disorder- induced large-scale inhomogeneities in density 
profiles persist if one relaxes the geometrical constraint by allowing the particles to move 
around defects? The purpose of our paper is to answer this question by considering two 
different realizations of bypassing in ASEP with disorder. In the first model we introduce 
long-range jumps weighted by a probability that decays with distance / as where 
0" > 1 to ensure finite current. The choice for this particular model is motivated by the 
fact that it may lack the phase coexistence in the presence of a single defect [25] . From 
the technical point of view, in the absence of disorder the model retains most of the 
characteristics of the short-range ASEP, but is more accurately described by the mean- 
field approach than in the short-range case [26]. The second model we consider is the 
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short-range ASEP in two dimensions, where disorder is present only in the direction of 
driving. This model is far from being just academic. It may be considered in a broader 
context of various transport phenomena in random media, ranging from a fluid flowing 
through a porous medium [271 [28], ^ traffic flow in the presence of obstacles [29] to, more 
recently, a transport of microfluidic droplets through various geometries [30] . Common 
to all these systems is the nontrivial interplay between driving, mutual interactions and 
the underlying geometry, which altogether may reduce the conductivity, in some cases 
even completely [27] . 

The remainder of the paper is organized as follows. Results for the long-ranged 
ASEP in one dimension are presented in section 2. The short-range ASEP in two 
dimensions is introduced in section 3, where we first analyse the current- density relation 
and then consider the anisotropic limit of hopping rates. A summary of results is 
presented in section 4. 

2. Bypassing by long-range hopping in one dimension 

We consider the model with N = pL particles distributed on L sites of a one-dimensional 
(ID)lattice with periodic boundary conditions, where each site is either occupied by a 
particle (r„ = 1), or is empty (r„ = 0). Dynamics is implemented by randomly choosing 
one particle during an infinitesimal interval dt and attempting to move it / sites to the 
right, where 1 < I < Lis taken from the probability distribution pi = / C/,_i(cr+l), 

being the partial sum of the Riemann zeta function. If the target site is empty, 
the move is accepted; otherwise, it is rejected. Additionally, we introduce A^d = cL (c 
being concentration) fixed but randomly distributed (defect) sites in a way that the rate 
of hopping is weighted by the additional factor r < 1 whenever particle jumps either 
to or from a defect site (see figure [1] for illustration). This particular choice of disorder 
fulfils two purposes: first, the particles can bypass the defect sites by the long-ranged 
jumps and second, the particle-hole symmetry is preserved. The latter is, in fact, not 
essential to our conclusions, but makes analysis of results easier. 





Figure 1. A schematic picture of the long-range hopping in disordered TASEP with 
periodic boundary conditions. Hopping rates are reduced by the factor r whenever the 
particle jumps from, or to a defect site (red colour). Note that not all the possible 
moves are sketched, but only those that illustrate the implementation of disorder. 
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2.1. Typical results of Monte Carlo simulations 

We are interested in the stationary density profiles (rj), i = 1,...,L, where (. . .) 
is defined as (...) = X]c(- • ^i.^) is the stationary solution of the 

corresponding master equation. Since the stationary solution is not known, the extensive 
Monte Carlo simulations were performed, representative results of which are presented 
in figures [2al and [2bl Figure [2a] shows the density profile typical of low densities p, where 
only the microscopic shocks appear. On the other hand, raising p above some threshold 
value pc induces a macroscopic shock (i.e. of size oc L), as presented in figure [2b] (due 
to the particle-hole symmetry, the same applies for lowering p below 1 — pc). 




i i 

Figure 2. Density profiles for (a) p ~ 0.1 and (b) p = 0.5 for a system of size L = 10^ 
witli r — 0.5, c — 0.5 and a — 1.8 obtained by Monte Carlo simulations {t = 10^ 
MCS/site) for a single disorder realization. 




<T.> 



Figure 3. (a) Density profile for a system of size L = 10"^ with r = 0.5, c = 0.5 and 
(7 = 1.2 obtained by Monte Carlo simulations [t — 10*" MCS/site) for a single disorder 
realization at density p — 1/2 and (b) the corresponding (normalized) histogram of 
densities (r^), i — 1 . . . ^ L which displays a single maximum around p — 1/2. 

So far, this is reminiscent of the behaviour found in the short-range case |23], except 
that Pc appears to depend strongly on a. This fact seems appealing as one may ask 
whether pc could reach 1/2 (resulting with absence of macroscopic shocks) by lowering 
0" > 1 (as it is the case for a single defect [25]). One choice of parameters for which 
this type of change may be observed in Monte Carlo simulations is presented in figure 
[3al along with the corresponding histogram, i.e. frequency distribution of densities (tj) 
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(figure [3bl) . This lack of phase coexistence may also be observed in the current- density 
relation j(p), where the usual plateau is missing (figure H]). 

On the other hand, finding pc for a general choice of parameters r, c and a proves 
to be an elusive task. Even in the short-range case one can only make an estimate, for 
example, using the so-called fully-segregated model [23] in which all the defect sites are 
in a consecutive order forming a large "bottleneck" . This is the approach we pursue in 
the following section. 
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Figure 4. Current-density relation (fundamental diagram) for various a {L = 10'', 
c = 0.5 and r = 0.5) obtained by Monte Carlo simulations (t — 10^ MCS/site) for a 
single disorder realization. For better comparison, the current has been scaled with 
the average hopping length Xl — i ■ Pi- 
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2.2. Mean-field approach 

In the fully-segregated model of Tripathy and Barma [23j, all the defects are in a 
consecutive order forming a large segment (say y) of size = cL, where the hopping 
rates are reduced by the factor r. Provided this domain is large enough, increasing of 
density p will eventually bring the stationary current j to its maximum value for this 
domain, j = r/4, which corresponds to density py = 1/2. Due to the conservation of 
current, j has to be matched with the current in the remaining segment (say x) 

Px(l-Px) =rpy(l -py). (1) 

Together with the conservation of particles, p^{l — c) + cpy = p, this gives the following 
(mean-field) estimate of pc 



Later, Krug [21] showed that ([2]) is actually the upper bound on the exact value of pc 
in the disordered model, while the lower bound is given by 



Pc > 7^ ■ (3) 
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Note that these two bounds meet in the hmit of infinitesimal concentration c — )■ (not 
to be confused with a single defect!) giving = {I — \/l — r) 12 7^ 1/2. 

For these results to be apphed to the original, fully disordered system, one argues 
that the largest "bottleneck" has diverging length / ~ InL in the limit L — )■ 00 [23], so 
that the current eventually establishes its asymptotic value joo = ''^Py(l ~ Py)- The same 
is true in the long-range case, except that the current joo gains an additional factor 
A = J2i ^ ' Pi^ Joo = A ■ rpy(l — py). Since this has no effect on ([1]), one expects the same 
Pc as in (I2]). This means that the absence of a macroscopic shock in figure l3al should 
be a finite-size effect and should be ruled out in the infinite system. However, as we 
shall see later, to observe it in Monte Carlo simulations would require enormous system 
sizes L (e.g. InL ^ 50, see I5b|) . [It is worth mentioning that such a slow approach to 
the thermodynamic limit was found in some other processes in driven diffusive systems 
[31j.] In what follows we therefore adopt the segregated model of Tripathy and Barma, 
but with / free parameter, i.e. not related to L. 




2000 4000 6000 8000 10000 "O 2000 4000 6000 8000 10000 

i i 



Figure 5. Density profiles in tlie segregated model for various I [L — 10"*, r — 0.5) 
obtained by Monte Carlo simulations [t = 10^ MCS/site) for a single disorder 
realization at density p = 1/2 for (a) a = 1.2 and (b) a = 1.05. In (a), dashed 
lines represent asymptotic values of pc (and 1 — p^) according to bounds ^ and (jS]) 
with c — l/L and / — 1000. In (b), onset of shock creation is captured for I that 
corresponds to system sizes unreachable in Monte Carlo simulations with full disorder. 

The results obtained by increasing / at two different values of cr, 1.2 and 1.05, are 
presented in figures [5a] and [5bl respectively. The density profiles are similar to those in 
the presence of a single defect, where shocks disappear with increasing r (short-range 
model) or with decreasing a (long-range model). The main difference is that in the 
single defect case, the system is able to attain the same value of the maximum current 
as in the pure model, which is done by building long-range correlations in the density 
profile. On the other hand, in the segregated or fully disordered models, this is not 
possible as the largest "bottleneck", whose length diverges with L, will force j to a 
value less than in the pure case, provided / is large enough. How large / is needed to 
establish the asymptotic regime is not an easy problem as it depends on the contribution 
to the current coming from the site-dependent corrections to px and Py (for example, in 
figure [5b| InL = 50 is still not enough to observe shock). Instead of calculating these 
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corrections explicitly, in what follows we estimate their leading contribution to the total 
current function of /. 

Let us enumerate lattice sites so that "bottleneck" occupies sites i = L — l + 1, . . . , L. 
In the long-range model, the current ji that satisfies local conservation law, 

= J.-l -J^, (4) 

is defined as the total current of all particles jumping from and over the site i, 

L-l 

J* = ^ ^ Pm+nin-mi^ — Ti+n))5l_^^i_^_n, (5) 
m=0 m+n<L 

where disorder is introduced through i equal to r or 1 depending on whether a pair 
of sites k, I contains at least one defect site or not, respectively. Note also that, due to 
the periodic boundary conditions, Tj±L = ^j, j = I, . . . , L. To calculate the current, 
we choose the site in the middle of the "bottleneck" (z = L — //2 + 1) and apply the 
mean-field approximation, (rjT,) ^ (rj)(rj). Assuming the following density profile 

(of, l<i< L-l 

gives four contributions to j (jxx, jxy, jyx and jyy) coming from the exchange of particles 
between x and y segments. Since we are interested only in how fast j converges with /, 
we can ignore p^'^ as they are bounded between and 1. It is then straightforward to 
show that 

Jxx <r(l-pi + 2-p2 + ... + (V2)PV2) -0(1) 

Jxy, Jyx < r(l ■ + 2 ■ pif2+2 + ... + {l/2)pi) ~ 0(/-('^-^)) (76) 

Jyy < 1 ■ pi+i + 2 ■ pi+2 + . . . + /p/ ~ 0(/~("~^)), (7c) 

so that contributions fl75|) and flTcj) decay very slowly when a is close to 1. 

To conclude, site-wise disorder in one dimension seems too restrictive, even if defects 
can be bypassed by the long-ranged jumps. As concluded by Tripathy and Barma, the 
reason for this is that disorder in one dimension generates "bottleneck" of slow sites that 
diverges in the limit L — )■ cxd. However, such geometric constraint is not to be expected 
in higher dimensions, which we consider in the rest of this paper. 



3. Bypassing by transverse hopping in two dimensions 

We consider = pL^Ly particles distributed on a two-dimensional (2D) quadratic 
lattice consisting of x Ly sites, where each site holds at most one particle at a time 
{Tij e {0, 1}). At any given moment, a randomly chosen particle attempts to move to 
the nearest site either in the direction of driving with the probability p\\ = p^ or in 
the perpendicular direction with the probability p± = Py so that p^ + 2py = 1. If the 
target site is empty, the move is accepted; otherwise, it is rejected. In other words, 
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the exclusion process is totally asymmetric in the x direction and symmetric in the ij 
direction with periodic boundary conditions assumed in both directions. To maximally 
ease bypassing, disorder is introduced by randomly choosing cL^Lj sites from where 
hopping in the direction of driving is forbidden, which could be imagined as if a "bond" 
between the two adjacent sites in x direction was broken. In other words, a particle at 
a defect site either moves transverse to the driving with probability Py (provided the 
target site is empty) or waits with the probability p^ until the next attempt. 

Two-dimensional exclusion process has been studied previously in various contexts. 
In [271 [2H], Ramaswamy and Barma have studied a partially asymmetric exclusion 
process with pup = Pright = w ■ {1 + g) and pieft = Pdown = w ■ {1 — g), where disorder 
is introduced by breaking "bonds" in both x and y directions with probability c. The 
underlying network created in that way is somewhat different from ours, because it 
creates backbends along which current flows against the driving field. Even closer to 
our problem is the work of Saegusa et al [29], who considered flow of particles in the 
multi-lane TASEP (Ly ^ L^) with fixed obstacles, where the usual plateau in the 
current-density relation was observed. Alexander and Lebowitz [32] considered the 
motion of a rod immersed in a fluid of interacting particles described by the symmetric 
simple exclusion process and observed a macroscopic region of low density behind the 
rod. Finally, large-scale inhomogeneities in two-dimensional driven diffusive systems, 
although not induced by disorder but akin to phase coexistence, have been investigated 
in many works, e.g. by Schmittmann et al [33]. 

3.1. Typical results of Monte Carlo simulations 

3.1.1. Current- density relation Our first task is to understand how the presence of 
disorder affects the stationary current in the x direction, which is defined as 

1 

Jx(p, Oi)/Ly = — ^p^{Tij{l - Ti+ij)) ■ ujij{a), i = l,...,L^, (8) 
y i=i 

where, for a given realization of disorder a, u}ij{a) = 1 if the bond that connects sites at 
(i, j) and {i + is present, and C0ij{a) = otherwise. Figure [6] shows j{p) obtained 
by Monte Carlo simulations for various concentrations c on 200 x 200 lattice with 
Px = 2py = 1/2. Compared to the pure case where jx{p)/Ly = Pxp{l — p), we find 
that, before reaching its maximum, the current can be quite well fitted to the parabolic 
shape, j{p) oc p(l — p), but the factor seems nontrivial and decreases with increasing 
concentration c. More importantly, one observes a plateau around p = 1/2, but its 
boundaries are less sharp compared to the one-dimensional model. 

A naive estimate of height of the plateau can be given by the average number of 
regular bonds per column times Px/4, (1— c)px/4, which is too high as seen from the third 
column in table [H This estimate can be improved by looking at the column with the 
smallest number of regular bonds, which can be either found directly in the particular 
realization of disorder or calculated using the extreme value theory [31] . Since we could 
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not find any reference in literature to point to tliis problem explicitly, we present details 
of this calculation in the Appendix and state here only the final result 



Ly — max < ^^[1 — i^ij(tt)] 
' U=i 

Ly - ai^ (c, Ly )7 - (c, Ly ) = cJ* , (9) 




0.08 



c = 20 % 
c = 30%. 
c = 40 % 




Figure 6. Current-density relation (fundamental diagram) for various concentrations 
c on 200 X 200 lattice with = 2py = 1/2, each obtained by Monte Carlo simulations 
(t = 10® MCS/site) for a single disorder realization. Dashed lines denote the best 
estimate of the maximum current given by the expression (jlOp . 



where 7 = 0.5772 ... is Euler-Mascheroni constant and ai^^c) and bL^{c) are given by 
the expressions flA.9p and flA.lOp . respectively, in the Appendix. [Note that uj* in ([9]) 
is independent of a since it is obtained by averaging over all disorder configurations 
with concentration c. However, since is generally large, this value is close to 00* {a) 
obtained by counting regular bonds explicitly from a particular realization of disorder 
a (as seen from the fourth and the fifth columns in table [I]).] Although the analytical 
expression ([9]) is closer to Monte Carlo data than (1 — cJp^/A, it poorly describes the 
maximum current overall, especially at higher concentrations c. The best estimate, as 
displayed in the sixth column of table [H is obtained if one recognizes that the largest 
contribution to the current will come from those sites (within a column) that have both 
inward and outward bonds, which gives 



max{j^{p,a)}/L 
p 




[Ly - aL,(2c 



Ly)]-pJ{AL 



(10) 



where 2 c 



1-fl 



is the probability of not having both inward and outward 



bond at the same site. It should be noted that the usefulness of this expression depends 
on the value of p^- in the strongly anisotropic limit where p^ ^ 2py, there will be an 
increasing contribution to the current coming from the particles entering the column at 
one site and then leaving it from another. Details of this analysis should be published 
elsewhere. 
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Table 1. Values of maximum current for various concentration c obtained from Monte 
Carlo simulations (Lx = Ly ~ 200, Px = 2px = 1/2, t ~ 10^ MCS/site) and compared 
to a naive guess (1 — c)px/4, the expression uj*{a)p-^/ (4Ly) obtained by counting bonds 
explicitly from a, the expression uj*Px/ (4iy) calculated from the extreme value theory 
and to the best estimate given by the expression pop . 



c 


Monte Carlo 


(1 - c)px/4 


u;*ia)-px/{ALy) 


ZJ* •Px/(4iy) 


expression (fTO|) 


0.1 


0.09255(5) 


0.1125 


0.10507 


0.10313 


0.09154 


0.2 


0.06953(7) 


0.1000 


0.09063 


0.09010 


0.06812 


0.3 


0.05089(2) 


0.0875 


0.07688 


0.07615 


0.04887 


0.4 


0.03553(5) 


0.0750 


0.06000 


0.06287 


0.03312 



3.1.2. Density profiles Having found some basic characteristics of the current- density 
relation for typical values of p^ and Py, we now turn to the investigation of stationary 
density profiles. Figure [7] shows the density profiles obtained by Monte Carlo simulations 
for various p, represented by the shades of blue ((Tjj) = 0), white ((Tjj) = p) and red 
{{r.;) = 1). 









■ \\:,- 



















■ 1 

1 ,':f' , 



1 

0.8 
0.6 
0.4 
0.2 




0-6|>!^iEKi'',™iitl,, , „ , , 



p = 0.4 





Figure 7. Density profiles for various p at c = 0.4 obtained by Monte Carlo simulations 
(Lx = Ly = 200, Px = 2px = 1/2, t — 10^ MCS/site) for a single disorder realization, 
represented by the shades of blue {{nj) — 0), white ((t.^) — p) and red {{rij) = 1). 



Inspecting profiles for increasing p, one clearly sees how local inhomogeneities 
present at lower p, start to grow until two large-scale phases occur for p around 0.4. 
This is even better seen from the corresponding histograms, which display either a 
single maximum at lower/higher densities (figure l8a|) or two maxima at intermediate 
densities (figure l8b|) . [It is also interesting to observe that even when p is small, density 
distribution is still not peaked around p] Whether this is a true transition or just a 
finite-size effect seems difficult to determine using only numerical approaches, as they 
are restricted to not too large system sizes. We therefore tackle the problem from a 
different angle by looking at the strongly anisotropic limit Px — 0, where a rather 
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Figure 8. Histogram of density profiles obtained by Monte Carlo simulations for 
a single disorder configuration at c = 0.4 (Lx = Ly = 200, = 2py = 1/2, 
t = 10^ MCS/site) displaying (a) single maximum at lower p and (b) two maxima 
at intermediate p. 



3.2. Mean- field approach in the limit p^ — t- 

In the limit 0, the jumps in the x direction are rare compared to the ±^ direction 
and Hi = Ylj '^ij niay be taken as a slow variable whose evolution in time is described 
by the following master equation: 

Ap(c, t) = Y^ [w{ci'+^ ^ c)P{ci'+\ t) - w{c ^ c!:'+')P(c, t)] , (ii) 

Lit/ 

t 

where C = {ni\i = 1, . . . , L^} and C^*'*'"'^ = {^i, . . . , rij ± 1, nj+i =F 1, . . . , til^}. The next 
question is how to properly describe transition rates W. Since jumps between columns 
are rare, we may assume that particles have enough time to distribute uniformly within 
each column, so that each site within a column i has a probability Ui/Ly of holding a 
particle. Under this assumption, the transition rates read 

W{C ^ C^^-') = p. . (^^c., j . ^ . (l - !|tl^ , (12) 

wicr' ^c)=p..(j: ■ ^ ■ - ^^^) ' (13) 

where ooi = ^jLi^ij- We have thus reduced the starting two-dimensional problem 
to the one-dimensional in which disorder is associated with each site through cuj. In 
this process, known as the misanthrope process (MP), the hopping rates depend on the 
position of the particular site as well as on the number of particles both at the departure 
site and at the target site. 

Let us write down the lattice equation for the average number of particles at site i 
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Note that this equation transforms into mean-field equation for the one- dimensional 
TASEP if we further assume (njrij+i) ^ and use (nj) — > {ni)/Ly = pi and 

uji — Ui/Lj = ri, which gives 



dp£ 
dt 



Px ■ ri_ipi_i(l - Pi) - ■ riPi{l - pi+i), 



where the are taken from the binomial distribution 



Pin = ri ■ L„ 



n 



(1 - c)"c^ 



(15) 



(16) 



with the mean and the variance given by {ri} = 1 — c and {rf} — {rj}^ = c(l — c)/Ly, 
respectively. 

To check how well the process defined by (fTT|) - (fT3ll describes the original data from 
the two-dimensional TASEP in the limit p^ — ?■ 0, we performed Monte Carlo simulations 
with the Ti derived from the same realization of disorder as in 2d TASEP and obtained 
stationary densities (ui). The results are presented in figures [9al and [9b] for p^ = 10~^ 
and Px = 10~^, respectively, together with the data from the corresponding 2d and Id 
TASEP. From these figures it is clear that the one-dimensional misanthrope process 
(MP) defined by f llip and (IT^ describes 2d data very well for several orders in p^, 
while the data from Id TASEP follow all the "peaks" but show generally less deviation 
from 1/2. As we raise Px further, it is evident that the picture of independent columns 
with uniformly distributed particles has to be abandoned. A crude estimate when 
this happens is when the relaxation time of fiuctuations within a column, r„ ~ L^, is 
comparable to the average time that a particle spends in one column, Tc ~ 1/Px- Since 
for the symmetric simple exclusion process z = 2, this gives p^ ~ L~^, which seems 
too small compared to the observed data from Monte Carlo simulations. The precise px 
when this happens is thus an open problem and further work in this direction is needed. 
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Figure 9. Column-averaged density profiles ((ri) = ^jZi{Tij) / Ly) obtained from 
Monte Carlo simulations of 2d TASEP for (a) = 10^" and (b) = 10"^ 
(Lx = Ly = 200, c = 0.2, p = 1/2, < = 10^ MCS/site) compared to density profiles 
obtained from the equivalent Id misanthrope proeess (MP) described by pT|) - (IT3)) and 
from the Id TASEP. 
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The "reduction" of the original two-dimensional problem to the one-dimensional 
allows us to invoke the same geometric argument as in Id when discussing phase 
coexistence. However, a striking difference between distribution (1161) of r^'s and the 
ones previously considered in the Id disordered TASEP ^231 |35] is in the variance. 
While in previous studies the variance always remained finite, here it scales as ^-/Ly 
meaning that in the limit Ly — )■ cxd it becomes improbable to observe an infinitely large 
domain of sites with, say, all rj less than some fixed value tq below average. To see this, 
let us recall the result obtained by the extreme value theory which gives the average 
length of the longest sequence of consecutive heads in coin tosses with P{head) = p 
[36], 



-f + \n[N{l-p)] 



1 

2' 



:i7) 



ln(l/p) 

where {. . .} denotes the average over all possible outcomes in N coin tosses. As 
mentioned earlier, this gives {/max} ~ Inlv when applied to the one-dimensional TASEP 
of size A^ with a binomial distribution of defects, where p is the concentration of defect 
sites. If instead of a binomial distribution one considers (fT6|) . the "bottleneck" may be 
defined as a domain of consecutive sites with all r^'s less then some fixed value ro < 1 — c, 
p = P{x < ro), as it was done by Krug in [21]. Approximating the distribution f|T6|) 
with the normal distribution with the mean and the variance given by /i = 1 — c and 
cr^ = c(l — c)/Ly, respectively, gives 



p = P{x < ro) 



1 + erf 



( 



ro- n 



V V2^ 

Now, if Lj is finite, both a and p are finite, which gives {/max} ~ Inl^x, i.e. an infinitely 
long "bottleneck" in the limit Lx — )■ oo. This result applies to the multi-lane TASEP 
with a finite number of lanes, for which this type of disorder should always induce phase 
separation. However, in the true 2d system of linear size L, where Ly L, the 

variance a — )■ in the limit Ly — )• oo. The error function erf(x) in ( ITSi) can be then 
expanded around oo. 



erf(a;) = 1 — 



2X' 



+ ... 



X — )■ OO 



(19) 



If, say, ro < /i = 1 



we may expand the error function around — oo using the fact 



that it is an odd function, erf (x) 
gives 



-erf(— x), which after some algebraic manipulation 



ln(l/p) = In 



2V^|ro- l + c| 



1 + c) 



(20) 



In the limit Ly — t- oo, the denominator in (|T71) diverges faster that numerator making 
{/max} vanish with increasing Ly. In other words, in the limit — it is as if 



I Similarly, one can set ro > jjL — 1 — c and then use the expansion (fT9| around oo. 
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"disorder averages itself" and the argument in favour of phase coexistence due to the 
large "bottleneck" fails. 

To check this in the original, 2d model, we present for illustration Monte Carlo 
simulations for the lattice size = Ly = 1000 with c = 0.4 (figure [TOj) . As it may be 



observed in figure llOal the current does not show the plateau, but instead can be rather 
well fitted to the expression i^{p,a) oc PxP(l — p) with the factor of proportionahty 
~ 0.5756 close to 1 — c = 0.6. In figure [TObl is displayed the corresponding histogram of 
density profile obtained for p = 1/2, which shows a single maximum with a Gaussian-like 
shape. 
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Figure 10. (a) Current-density diagram obtained by Monte Carlo simulations of 2d 
TASEP for a single disorder configuration at c = 0.4 (px = 10^^, ~ Ly ~ 1000 and 
t = 10^ MCS/site). Solid line is the best fit to the expression j{p) — const, x p(l — p) 
with const, w 0.5756 close to 1 — c = 0.6. (b) The corresponding histogram of density 
profile for (0 = 1/2. 



4. Conclusion 

In this work we studied the possibility of avoiding macroscopic phase separation in 
the exclusion processes with site-wise quenched disorder by allowing particles to bypass 
defect sites. Bypassing was implemented in a natural way, either by extending the range 
of hopping or by increasing the dimensionality. In the first example, the performed 
Monte Carlo simulations indicated the absence of the macroscopic shock provided the 
hopping length / is taken from the probability distribution obeying the power-law 
Pi ~ and a is close to 1. This was a rather surprising result, since the origin 

of shocks in one dimension is related to the fact that disorder creates a cluster of slow 
sites whose length diverges in the thermodynamic limit. However, by studying the 
entirely segregated model in which all defects are in the consecutive order forming a 
large "bottleneck", we showed that the absence of shocks is a finite-size effect which 
decays very slowly with the size I of the "bottleneck" and thus persists for enormously 
large system sizes, e.g. / ~ InL ^ 50. 

In the second example, we considered the 2d exclusion process, which is totally 
asymmetric in the longitudinal direction and symmetric in the transverse one, with 
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disorder introduced through obstacles in the driving direction. We have shown that 
a rather simple quantitative description of the process leads to the absence of phase 
separation in the strongly anisotropic limit, where the hopping rate in the direction of 
driving is much smaller than in the perpendicular direction, px ^ 2py. In this limit, the 
original two-dimensional exclusion process reduces to the one-dimensional misanthrope 
process, in which disorder enters only through the fraction of non-defect sites present 
in each column of the original 2d TASEP. By the central limit theorem, however, the 
probability distribution of this fraction has a variance which decays as so that 

in the limit Lx ~ — )■ oo the otherwise diverging size of the largest "bottleneck" 
vanishes resulting in the absence of a macroscopic phase separation. On the other hand, 
if Lj remains finite while Lx — ?■ oo, as in the multi-lane TASEP with a finite number of 
lanes, the argument in favour of bottleneck-induced phase separation should be revoked. 
Obviously, the same argument remains valid for more general types of disorder, provided 
it obstructs flow only in the direction of driving. 

Although the general properties of 2d disordered exclusion process remained out 
of our scope, some interesting open questions emerged. First is the upper value of Px 
for which the above mapping applies. Second is the mechanism of phase separation for 
Px outside this limit and related to that is the finding of threshold density pc- Last, 
it would be interesting to examine the case where the bond disorder, whose realization 
may be found in some microfluidic geometries [30], is replaced by the site disorder, 
appropriate for modelling flow of a fluid through porous media or pedestrians avoiding 
fixed obstacles. 
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Appendix 

If Xi denotes the total number of broken bonds in the i-th column, Xj = ^^Ii(l —oJij), 
then the corresponding probability distribution is the binomial distribution 



where c is the probability of finding a broken bond at site {i,j)- Let F{x) = P{x < X) 
be the corresponding cumulative distribution and x* its right endpoint, x* = sup{x : 
F{x) < 1}. We are interested in obtaining the maximum value of {Xi, . . . , Xl^} as 
Lx — )■ oo. 



where -^^ means convergence in probability, since P{max{Xi, . . . ,Xm} < x) = F^{x) 
is degenerate in the limit m — )■ oo as it converges either to for x < x* or 




(A.l) 



max{Xi, ... , Xl^} x*, 



(A.2) 
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to 1 for X > X* . We therefore seek a sequence of positive and real such 
that lirrim-^ooF'^ {cimX + &m) = G{x) exists, where G{x) is called the extreme value 
distribution. A sufficient condition for that is von Mises' condition [34], which states 
that if F"(x) exists and F'(x) > for x < x* then 



Urn F^{amX + hm) = exp 



(1 +7'x)-^/^'l , l + 7'x>0. 



where 7' is given by 



lim 



[1 - F{t)]F"{t) 



-7 



1. 



(A.3) 



(A.4) 



[F'{t)Y 

Moreover, hm = U{m) and = m,U'{m), where U{m) is the inverse function of 
1/(1 — F{x)). For 7' = 0, the right-hand side of ( 1A.3I) should read exp(— e^"^). To 
apply this condition, we approximate a binomial distribution with the normal N{fi, o"^), 
where fi = cLy and o"^ = Lyc{l — c). Numerical error in doing so need not worry us 



since for cumulative distribution it is of the order of 1/ \/Ly (Berry-Essen theorem) and 
Ly is large. Then it is an easy exercise to show that the cumulative distribution F{x) 
and its inverse U (x) are given by, respectively, 

1 



Fix) 



U{x) 



1 + erf 



fi + V2a2 ■ erf 1 - - 

X 



X > 1, 



(A.5) 



(A.6) 



where erf(x) and erf~^(x) are the error function and its inverse, respectively. Inserting 
F{x) and its derivatives in ( ]A.4p we obtain 7' = 0, i.e. G{x) = exp(— e~^) (Gumbel 
distribution). The mean and the variance of the Gumbel distribution are given by Euler- 
Mascheroni constant 7 = 0.5772. . . and 7r^/6, respectively, which gives the mean and 
the variance of x* 



*\2 



6 



where and hi^ are given by 
^lAc, Ly) 



c)U 



exp 



erf'^ 



1-A 



6i,(c, Ly) = cLy + ^2c{l-c)Ly ■ erf-1 ' ' 

For > 1, erf-^(l - can be expanded around Lx — ?■ 00, which gives 

4c(l - c)Ly 



«Lx(C, Ly 



27rlnL^ — 7rln2 



TT 



{C,Ly) ^ 



cLy + < c(l — c)Ly 



1/2 



In 



2tt 



— Inln 



2n 



1/2 



(A.7) 
(A.8) 

(A.9) 
(A.IO) 

(A.ll) 
(A.12) 
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